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The classical theory of Riemann ellipsoids is formulated naturally as a gauge 
theory based on a principal G-bundle V . The structure group G = SO (3) is the 
vorticity group, and the bundle V = GZ/_|_(3, R) is the connected component 
of the general linear group. The base manifold is the space of positive-definite 
real 3x3 symmetric matrices, identified geometrically with the space of inertia 
ellipsoids. The bundle V is also a Riemannian manifold whose metric is deter- 
mined by the kinetic energy. Nonholonomic constraints determine connections 
on the bundle. In particular, the trivial connection corresponds to rigid body 
motion, the natural Riemannian connection - to irrotational flow, and the in- 
variant connection - to the falling cat. The curvature form determines the 
fluid's field tensor which is an analogue of the familiar Faraday tensor. 



1. INTRODUCTION 

A Riemann ellipsoid is a constant density fluid bounded by an ellipsoidal 
surface and constrained to a linear velocity field. In addition to rotational 
and vibrational motion, Riemann ellipsoid dynamics admits vorticity de- 
grees of freedom. The vorticity allows the moment of inertia to be varied 
continuously between its rigid body and its irrotational fluid limits. It is 
this flexibility that enables the theory to model diverse rotating systems 
[1, 2, 3, 4, 5]. 

An ellipsoid is defined by a positive-definite real symmetric 3x3 matrix q; 
the points X — [X\, X%, A3) of its surface in Euclidean space are solutions 
to the quadratic equation 



Qu-XiXj = 1- (1) 



In particular, when q = diag(a^, a|, &§), the surface is an ellipsoid with 
major axis half-lengths equal to a\, et2, 0,3 and the ellipsoid's principal axes 
are aligned with the Cartesian axes. The space of all ellipsoidal surfaces is 
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identified with the six dimensional manifold of such matrices q: 

Q = {qeM 3 (R) | q t =q,q>0}. (2) 

The configuration space for Ricmann ellipsoids is the general linear group 

V = {£ G M 3 (R) | dct£ > 0} = GL+(3,R), (3) 

because only linear motions are allowed. The elegance and beauty of the 
theory of Ricmann ellipsoids are derived ultimately from the group struc- 
ture on the configuration space V. 

The general linear group transforms an ellipsoid into another ellipsoid 
that is deformed and rotated. If q G Q and £ G V, then this spatial 
transformation corresponds to q i — ► In particular, the unit sphere 

with q = I is changed into the ellipsoid with q — Thus, the natural 
mapping tt from the group V onto the ellipsoidal space Q 

tt-.v — > q 

£ — > = (4) 

is surjective (onto). 

The group V also acts on itself. The group elements can be multiplied in 
two ways, on the left or on the right. The left and right actions are distinct 
when the group is not Abelian. They define two different transformations 
of the group on itself. Left multiplication of elements £ G V by r € 50(3) 

^ LrZ = r£ (5) 

corresponds to physical rotations r in three-dimensional Euclidean space. 
In contrast, right multiplication of £ G V by g g G = 50(3), 

£~R g t = £g- 1 , (6) 

corresponds to vortex motion. Although the groups G and 50(3) are math- 
ematically isomorphic, their distinct actions on V imply different physical 
interpretations. 

With respect to left multiplication by the rotation group element r G 
50(3), the ellipsoid with inertia tensor q = = is transformed 

into the rotated ellipsoid with inertia tensor 7r(Lr£) = r££*r* = rqr t . But 
right multiplication by the vorticity group clement g G G leaves the inertia 
ellipsoid invariant, ir(R g ^) = £.g _1 g£ t = q, since g^ 1 = g*. This invariance 
is expressed more elegantly by the composition of mappings 



tt o R g = 7T, for all g G G. 



(7) 
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Hence G is a group of motions internal to the ellipsoidal surface. 

Because of the right invariance of the surjective mapping n with respect 
to the group G, the configuration space V is a principal fiber bundle over 
the base manifold Q with structure group G [6, 7]. It is the purpose of 
this article to show that the classical theory of Riemann ellipsoids is ex- 
pressed naturally in the mathematical framework of the bundle V. The 
bundle V has two additional structures that are important physically. V 
is a Riemannian manifold whose metric is determined by the kinetic en- 
ergy in three-dimensional Euclidean space. V may also be equipped with 
a differential structure, or connection. The nonholonomic constraints to 
irrotational flow and the falling cat problem correspond to the Riemannian 
connection and the invariant connection, respectively. 

Recently Littlejohn and Reinsch [8] and earlier Shapere and Wilczek [9] 
presented a gauge theory of rotating systems in which the gauge group 
is the rotation group, and the gauge conserved quantity is the angular 
momentum. This article describes a fundamentally different model in which 
the gauge group is the vorticity group and the gauge conserved quantity is 
the Kelvin circulation (the angular momentum is also conserved). 

The mathematical theory of bundles is the setting for the gauge theo- 
ries of fundamental forces [10, 11]. For the Riemann ellipsoid model, the 
gauge group is the vorticity group G — SO (3), and, therefore, it is a non- 
Abelian gauge theory. The situation is reminiscent of the Weinberg-Salam 
theory of the electroweak interaction [12] for which the structure group 
is SU(2) x U(l). An important mathematical difference is that the base 
manifold for the electroweak problem is Minkowski space-time, while the 
Riemann ellipsoid base manifold is the space of ellipsoidal surfaces. Al- 
though Cartesian coordinates (i, x, y, z) are convenient for space-time, a 
coordinate system for Q requires Euler angles to describe the ellipsoid's 
orientation. This technical complexity can be overcome by using Cartan's 
method of moving frames [7]. For the rotation group SO (3), the angular 
momenta serve as the moving frame, and, for the gauge group G, the Kelvin 
circulation is used. In this sense, the differential geometry of the Riemann 
ellipsoid theory is more complicated mathematically than the space-time 
based theories of fundamental physics. 

In §2 the theory of rotating triaxial rigid bodies is presented using group 
theoretic and geometrical methods that generalize to the Riemann ellip- 
soid problem. The rigid body configuration space is SO(3) that is both 
a Lie group and a Riemannian manifold. Applying the Cartan method of 
moving frames, the velocity of a rotating body is expressed in terms of the 
angular momentum which is an invariant vector field. The kinetic energy 
determines a metric on the rotation group. The Euler rigid body equations 
are obtained from Lagrange's equations. 
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The Ricmann ellipsoid theory is presented in §3. The interrelated Lie 
group, bundle, and Riemannian structures of the configuration space V 
are reviewed and the equations of motion are derived in the Lagrange for- 
malism. §4 is the heart of the paper. Here the differential geometry of 
bundle connections is used to show that the nonholonomic constraints to 
irrotational flow and the falling cat are realized by mathematically natural 
connections. 

2. RIGID BODY ROTATION 

Although the rigid body problem is solved in basic mechanics textbooks, 
the method of solution does not generalize easily to higher dimensional 
non- Euclidean configuration manifolds, e.g., the Riemann ellipsoid space V. 
The source of the limitation in the conventional technique is the use of Euler 
angle coordinates for the configuration manifold 5*0(3). In this section the 
rigid body problem is solved using the method of moving frames. Cartan's 
more powerful technique enables the solution of mechanics problems in 
higher dimensional non-Euclidean spaces. 

In the conventional solution of the rigid body problem, Euler angles 
provide coordinates for the configuration space and a basis for each tangent 
space. For the rigid body problem, the natural physical and geometrical 
choice for a tangent space basis consists of the angular momentum vectors. 
No coordinate system can define directly the angular momenta because 
their Lie brackets do not vanish. Cartan's moving frames provide each 
tangent space with bases that vary from point to point and cannot be 
expressed consistently in a fixed coordinate chart. 

2.1. Tangent Vectors 

A vector Vr tangent to the group manifold 5*0(3) at the point R de- 
termines a unique directional derivative as follows: if / is a real-valued 
differentiable function defined on an open neighborhood of the point R in 
the rotation group, then the real number {Vnf){R) is defined as the deriva- 
tive of / in the direction Vr evaluated at R G 50(3). A vector field V on 
50(3) is defined by a smooth assignment of a directional derivative Vr to 
every point R G 5*0(3). Hence a vector field V may be regarded as a first 
order differential operator on 50(3): if / is a smooth function on 50(3), 
then Vf is also a smooth function defined by (V f)(R) = (Vkf)(R), for 
all R G 50(3). The space of all vector fields on 50(3) is identified with 
the linear derivations, i.e., the space of all mappings V from the space of 
smooth functions on 50(3) into itself that satisfy the linearity and Leibniz 
conditions, 



V(af + bg) 



= a(Vf) + b(Vg) 
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V(fg) 



f(Vg) + (Vf)g, 



(8) 



where / and g are smooth functions on 5*0(3) and a, b are real numbers. 

The Lie bracket of two vector fields is the commutator of their corre- 
sponding derivations. Since it is a first order differential operator, the 
bracket of two vector fields is a vector field. In this way, the space of all 
vector fields is an infinite-dimensional Lie algebra. The components of a 
vector field depend on the choice of a coordinate system for 50(3). If Eulcr 
angles (fl 1 ,^ 2 ,^ 3 ) are taken, then the components of a vector field V are 
given by applying V to the coordinate functions, V{9 1 ) — V 1 . Expanded 
into components, the vector field is 



because 30 /d6 i = 5{. 

There are two principal advantages to regarding vector fields as first order 
differential operators. The first is that the Lie bracket may be defined 
easily. The second is that the differential operator V is independent of 
the coordinate system, although its components are coordinate dependent. 
If a second coordinate system is chosen for 50(3), say (ip 1 , ip 2 , V' 3 ), then 
the vector field components change to V{tp l ) — v % . But the differential 
operator, expressed in the new coordinate chart, stays the same, 



since the components of a vector transform as the Jacobian, v l (dO? / dip 1 ) = 
VK 

It is a tactical error to work directly with Euler angles or any other 
non-Euclidean coordinate system. Even if results are obtained successfully, 
the derivations are unilluminating, and can make trivial results seem very 
complicated, c.f., Goldstein's derivation of the Euler rigid body equations 
[13]. Instead it is best to use consistently a Euclidean coordinate system. 
Consider the nine-dimensional space M^M) of all 3 x 3 real matrices. The 
rotation group 50(3) is viewed as a closed, bounded three-dimensional 
surface in M3(R). The nine entries of each matrix provide a convenient 
global coordinate system for M3(R). In terms of this coordinate system, 
the points of the 50(3) surface satisfy two polynomial conditions that are 
naturally expressed in matrix notation 



There are six independent quadratic polynomial conditions on the 9 matrix 
coordinates Rij, viz., J2k RkiRkj = for i < j. Hence the surface 50(3) 




(9) 




(10) 



50(3) = {R e M 3 (M) | R*R = I, dcti? = 1} . 



(11) 
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is 9 — 6 = 3 dimensional. When a rigid body is rotating, its motion defines 
a smooth curve t i — ► R(t) in the surface of 50(3). The tangent or velocity 
vector to this curve in the surface 50(3) is given in the matrix coordinate 
chart of M 3 (R) by 

V {t) = Y J R l3 {t)^-. (12) 



2.2. Left and Right Invariant Vector Fields 

Among all possible vector fields, the invariant ones are singled out for 
their mathematical simplicity and physical relevance. In the following dis- 
cussion, the left invariant vector field at a point R(t) is shown to correspond 
to the angular velocity in the laboratory (inertial) frame, and the right in- 
variant vector field at a point R{t) is shown to correspond to the angular 
velocity in the body-fixed (intrinsic) frame. 

For each element of the Lie algebra of the rotation group, a left invariant 
vector field and a right invariant vector field can be defined. Let 50 (3) 
denote the Lie algebra of the rotation group 190(3), 

so(3) - {lie M 3 (R) I n* = -Q}. (13) 

As it is customary, the upper case SO (3) denotes the group while the lower 
case 50 (3) denotes the algebra. The rotation of a rigid body t i — ► R(t) 
naturally defines a Lie algebra clement at each point. Since R t R = I 
at each time t, its vanishing derivative with respect to t implies that the 
matrix 

n L (t) = R l R (14) 

is antisymmetric. Vl L has the physical interpretation of the angular velocity 
in the laboratory (inertial) frame. In this article superscript L will indicate 
a quantity in the laboratory (inertial) frame. 

A basis for the three dimensional vector space SO (3) is the set of an- 
tisymmetric matrices ti for i = 1,2,3, with matrix elements given by 
(ti)jk = £ijk- If ^ L £ 50(3) is expanded in this basis, il L = 
u> L is interpreted as the laboratory frame angular velocity vector. 

When Q £ 50(3) is independent of time, define (£n) e to be the tangent 
to the curve 9 i — ► exp(6>f2) in the group 50(3) that passes through the 
group identity e when — 0. For example a rotation about the 3-axis is 
given by 

/ cos 9 sin 9 \ 
exp(<9e 3 ) = -sin0 cos9 . (15) 
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A left invariant vector field is defined by spreading this tangent vector 
around the group manifold. At the point R 6 5*0(3), define (£oJ.r to be 
the tangent to the curve 9 i — > R ■ cxp(6>f2) that passes through R G 5*0(3). 
Similarly a right invariant vector field (3?n)fl is the tangent to the curve 
9 i — > cxp(— 9Q) ■ R, see Figure f . 

As first order differential operators on 50(3), the invariant vector fields 
are given by 



(H a ) R - (ffl) y ^- (Km) 



(K n ) R = -(fiflJyJL, (16b) 

for each £1 e SO (3). These derivations are called Lie derivatives. 

The two vector fields are tangent to the three dimensional surface of or- 
thogonal unit determinant matrices contained within the nine dimensional 
manifold of all 3 x 3 real matrices. The velocity vector of a curve 1 1 — > R(t) 
may be expressed at each instant in time as the value of left invariant vector 
field at the point R(t) 

v(t) = E air = ? (i2fiL(t)) «alr = {Lql u )r w ■ w 

As Q ranges over the basis ti of SO (3), the left and right invariant vector 
fields provide separate bases for each tangent space to the group manifold 
50(3). To see the relation between the two invariant bases, note that the 
curve whose tangent at R is a left invariant vector field may be expressed 
alternatively as 

9 i — ► R-exp(9n) 
= Rcxp(9fl)R- 1 R 

= exp(6Ad R Cl)R, (18) 
which is a curve tangent to a right invariant vector field: 

(£>n)R = -(%Ad R n)R. (19) 

Here the adjoint transformation Ad is introduced: For each R e 50(3) 
and il G SO (3), define the Lie algebra element Ad/jfi = RflR -1 in SO (3). 
The adjoint mapping U i — > AduCt is equivalent to ordinary rotation of the 
corresponding velocity vector lj i — ► Rw. The equivalence depends on the 
determinant of a special orthogonal matrix being equal to +1, and, hence, 
u) is an axial vector. 
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If fl L (t) is the laboratory frame angular velocity, fl(t) — Ad R ^ fl L (t) = 
RR* is the angular velocity in the rotating frame. Hence, the velocity in 
the laboratory frame is expressed in the rotating frame by 

V(t) = {&n L (t))R(t) = -(^o(t))fl(t)- (20) 

When ft = ^i^iti is expanded in the usual basis for S0(3), the vector 
u) = Riv L is called the angular velocity vector in the rotating frame. 
The left and right invariant vector fields are special for various reasons: 

1. Each is a representation of the Lie algebra, i.e., 

£[Qi,n 2 ] = [£ni,£n 2 ] ( 21a ) 
Vi,fi2] = PW^nJ. (21b) 

The invariant vector fields are finite-dimensional subalgebras of the alge- 
bra of all vector fields. Both the left and the right invariant vector fields 
are faithful, but reducible, representation of the Lie algebra SO (3). With 
respect to the Haar measure, each Lq and is a skew-adjoint operator 
on the space of squarc-intcgrablc functions on 5*0(3). In particular, the 
Hermitian operators 

Ij = ify (22) 

are the quantum mechanical angular momenta in the laboratory (inertial) 
frame and the body-fixed (intrinsic) frame, respectively. 

2. As Lie derivatives, they are invariant with respect to the group mul- 
tiplication. For each g in the SO(3), let L g and R g denote left and right 
multiplication, 



L g {R) = gR (23a) 
R g (R) = Rg-\ (23b) 

for all R in 50(3). Both multiplications are group homomorphisms, e.g., 
L g Lh = L g h for all <?, h in the group. The group homomorphism L g maps 
a neighborhood of R £ SO (3) to a neighborhood of gR <G 5*0(3). The 
curve 9 i — ► Rcxp(8il) through R is mapped by L g to the curve i — ► 
gRexp(8Q) through gR. Hence the left group multiplication induces a 
mapping, denoted by (£ s )*, from the tangent space at R to the tangent 
space at L g (R) = gR. But the tangent vector (L 9 )*(£q)k is none other 
that (£jn) g R. In this sense the left invariant vector field is invariant with 
respect to the left group multiplication. Similarly, a right invariant vector 
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field is invariant with respect to right group multiplication 



{L g )*(Ln) R = {&n) g R (24a) 
(R g )*(Rn)R = {$a)Rg-i. (24b) 

Unless a rigid body is rotating with a constant angular velocity tt L , the 
velocity V(t) will not be left invariant along the trajectory t i — > R(t). 



Although at each instant in time the velocity vector along a trajectory 
may be expressed as the left invariant vector £^1^, it keeps changing 
with time. A similar remark applies to the right invariant vector 3lji(t) . 

3. What happens when (R g )* is applied to a left invariant vector field? In 
this case, the curve 6 i — ► Rcxp(8il) through R is mapped by R g to the 
curve 9 i — > Rexp(8fl)g^ 1 through Rg^ 1 . But 

flexp^fi)^ 1 = R(g- 1 g)exp{6{l)g- 1 
= {Rg-^cMOgVg- 1 ) 
= {Rg- l )exp{6AA g Sl). (25) 

Therefore (R g )* maps the left invariant vector field corresponding to fl 
at R to the left invariant vector field corresponding to Ad s (Q) at Rg~ x , 
see Figure 2. A similar argument holds for the right invariant vector fields 
acted on by (L g )* 

(R g )*(Ln)R = (^Ad 9 n)ij s -i (26a) 
{LgURnh = (^A dg n) 9 R- (26b) 

In particular, if g = R(t) is the actual rotation of a rigid body, then at 
each instant 

(R R{ t))*V(t) = 0Gn(t))e. (27) 

4. The set of left invariant vector fields {Lgj , £g 2 , Lg 3 } provides a basis 
for each tangent space. It is called a "moving frame." The velocity is 
V(t) = Y^^i (t) (^ejit(t)- Similarly the set of right invariant vector fields 
{Dlg 1 ,3le 2 ,3le 3 } is another moving frame and V(t) = — (^ejfl(t)- 

2.3. Riemannian Structure 

The group manifold 5*0(3) is a Riemannian manifold, i.e., there is a 
positive-definite metric defined on it. The metric depends on the moments 
of inertia of the rigid body and is chosen so that the rigid body's kinetic 
energy is proportional to the squared length of the velocity. Thus the 
metric on 5*0(3) is inherited from three-dimensional Euclidean space. 

Let X Q denote the position vector in an inertial frame for particle a of 
mass m a in a system of many particles (a denotes the particle index) . Let 
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M = ^2 m a denote the total mass of the system and let Rq be a convenient 
unit of length. The dimensionless quadrupolc-monopole tensor q(X) for a 
collection of particles is defined by 

(MR 2 /5) q{X)ii = Y, m a X ai X aj , (28) 

a 

where the sum is over the particle index a and X a i, i = 1, 2, 3, denotes the 
Cartesian components of the position vector. For a continuous fluid the 
sum is replaced by an integral over the mass density distribution. 

For rigid rotation there is a curve t i — ► R(t) in the rotation group 
for which i?(t)X a (t) = x Q is independent of the time t for all particles 
a. The motion is a collective rotation because all the particles move in 
an orchestrated manner according to the same linear transformation R(t). 
As a real symmetric positive-definite matrix, q(X) may be diagonalized 
by a rotation and its eigenvalues are real positive numbers. Hence, the 
reference distribution x Q may be chosen so that its quadrupole-monopole 
tensor q(x) = Rq(X)R is a diagonal matrix, 

(MRl/5) q(x)ij = Y m a x ai x aj = {MRl/5)a% (29) 

a 

or q(x) = A 2 = diag(a^, a|, a|). For a uniform density ellipsoid, a* is the 
half-length, in units of Ro, of the i th principal axis. 

The velocity vector for each particle is U Q = i?'x Q . The kinetic energy 
of the system is 

t = £^u Q .u Q 

a 

= (MR 2 /10) Y^'kil'ki 

ij 

= -{MRl/lQ) tr(A 2 n 2 ) 

= ^{huif + h^ + hojl), (30) 

where the moments of inertia of a rigid ellipsoid are I\ — (MR 2 ) /b){a2+a\), 
etc. The last line is proven using the identity (fi 2 )ij = UiUj — 5ijU ■ u). 
At each point R G SO (3), define the metric g as the bilinear form 

gRKOnJfl, (XnM - -tr(0!A 2 Q 2 ) (31) 

for each pair of tangent vectors (^njfi, (Rq 2 )r at R where Oi, O2 G SO (3). 
This metric is manifestly right invariant which is the expression of its rota- 
tional invariance in the laboratory frame. Moreover, the rigid body kinetic 
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energy is proportional to the squared length of its velocity 

T=(MR 2 /10)s m (V(t),V(t)). (32) 

2.4. Lagrange's Equations 

The dynamical equations for the rigid body may be derived from its La- 
grangian. The Lagrangian is just the kinetic energy and the action is its 
integration along paths between two fixed points in 5*0(3). Rigid body 
motion corresponds to the path that minimizes the action with respect 
to all such paths. If Eulcr angles are introduced, the calculation of the 
minimization equations is tedious. It is better to enlarge the space from 
5*0(3) to the 12-dimensional Euclidean configuration space whose coordi- 
nates are the nine matrix entries of R and the three body-fixed angular 
velocity components of ft [14]. Since the allowed paths in the enlarged 
space must correspond to paths in the original configuration space 50(3), 
they must satisfy the nine constraint equations 

= RtR-I 

= RR 1 — ft. (33) 

Introduce nine Lagrange multipliers in the form of real 3x3 symmet- 
ric a and antisymmetric /i matrices. Then the Lagrangian on the twelve 
dimensional configuration space is given by 

L[R, R, 0, fi] = T(Sl) + ±tr[cr • (R f R - I)} + ±tr[/i • (A/Z* - fi)]. (34) 

This is just a quadratic function of the 12 coordinates and their deriva- 
tives. To find the equations of motion, note first that the Lagrangian is 
independent of f2. Therefore, the Euler-Lagrange equation simplifies to 
dL/duJi = 0, or 

dT 

Wj = ~ £ijk d^~ k (35) 

The derivative of the kinetic energy with respect to the angular velocity is 
the angular momentum in the body-fixed frame 

L k = u>kh ■ (36) 

This is an axial vector and is naturally represented by an antisymmetric 
matrix L = J2L k c k ; therefore, the Lagrange equation for the ignorable 
coordinates says that \i is the negative of the angular momentum 



jJL = —L. 



(37) 
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Next, the Lagrange equations for the Rij give the matrix equation 

(i- R + 2^- R = R- (a + cr*) (38) 

On the constraint hypersurface, the matrix 

(i + 2n-n = R-(a + a t )-R t (39) 

is symmetric, or 

A = [n, M ]. (40) 

Therefore, the dynamical equations for a rigid body are a finite dimensional 
Lax system [15, 16, 17] 

j t L=[n,L}. (41) 

Expressed in vector form, this becomes the usual Euler equation for the 
precession of the angular momentum vector 

^ = -<jJ x L. (42) 
dt K ' 

In addition to rotational symmetry, there is a symmetry of the rigid 
body due to invariance with respect to certain rotations in the body-fixed 
frame. While rotational symmetry follows from the invariance of the metric 
with respect to right multiplication, a body-fixed symmetry corresponds to 
metric invariance with respect to left multiplication. For g G 50(3), left 
multiplication by g of the metric yields 

= -tr (Ad g n 1 A 2 Ad g Q, 2 ) 

= -tr(fii {g- 1 A 2 g) Q 2 ) . (43) 

When g and A 2 commute, A 2 g = gA 2 , the metric is invariant with respect 
to left multiplication by g. Let H denote the subgroup of 5*0(3) consisting 
of the rotations g that commute with A 2 . If a curve t i — > R(t) minimizes 
the action, then, for each g G H, so must the curve t i — ► gR(t) that 
has the same action. Hence, the dynamical system factors through the 
projection onto the right coset space SO(3)/H, and the solutions are curves 
t i — ► HR(t) in this coset space. Since the metric is left iJ-invariant on 
50(3), it induces a well-defined metric on the right coset space SO(3)/H. 
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There are three possibilities for the subgroup H . For a triaxial rotor, the 
generic situation, the three principal moments are different, say a\ > a 2 > 
a 3 > 0, and the subgroup H is the dihedral group 

D 2 = {m = diag(mi, m 2 , m 3 ) | m t = ±l,det(m) = 1}. (44) 

It is a discrete group with four elements corresponding to rotations by ir 
about each principal axis. When the rotor is axially symmetric, two mo- 
ments are equal. For prolate (a\ > a 2 = a 3 > 0) and oblate (ai — a 2 > 
a 3 > 0) shapes, the subgroup H = 0(2), and it consists of all continuous 
rotations about the symmetry axis plus a rotation by it about a perpendic- 
ular nonsymmetry principal axis. In the case of a sphere (a\ = a 2 = a 3 ), 
H = SO (3) and the coset space is a fixed point. 

The coset space is geometrically interpreted as the space of orientations 
of an ellipsoidal body. Consider a rigidly rotating triaxial ellipsoid with 
fixed axis lengths oi > a 2 > a 3 > 0. Let O denote the closed submanifold 
of Q that corresponds to the orientations of this ellipsoid, 

O = {q = R t A 2 ReQ,Re 50(3)} . (45) 

Although each rotation R defines a point of O, the correspondence is not 
one-to-one. If R\ and R 2 define the same oriented ellipsoid, q = R\A 2 Ri = 
R\A 2 R 2 , then R 2 R\ commutes with the diagonal matrix A 2 , i.e., R 2 R\ S 
D 2 and these two rotation group elements are in the same right coset, 
D 2 R 2 = D 2 R 1 . Thus the space O of orientations of a triaxial ellipsoid is 
canonically diffcomorphic to the right coset space 50(3)/ D 2 . Moreover, 
a rigid rotation of the ellipsoid, q i — > r l qr for r G 50(3), is equivalent 
to right multiplication in the right coset space, D 2 R i — > D 2 Rr. If the 
ellipsoid is axially symmetric, then the configuration manifold is the right 
coset space 50(3)/0(2). For a sphere, the configuration space is a single 
point. 

2.5. Vibrations 

If the three axis lengths vibrate, then the number of degrees of free- 
dom must be increased at least by three. To model rotational and vibra- 
tional motion, a minimal extension of the theory is attained by requiring 
that there exists a curve t i — > R(t) in the rotation group and a curve 
t i — > A(t) in the space A of positive-definite diagonal matrices for which 
R(t)~K a (t) = x a (i) depends on time t, but A(t)^ 1 x a (t) = y a is indepen- 
dent of time. Since all the particles move together according to the same 
linear transformations R(t) and A(t), the motion is a collective rotation- 
vibration. As for rigid rotation, the orthogonal matrix R(t) is chosen so 
that q(x(t)) = A(t) 2 is a diagonal matrix. The diagonal matrix A(t) makes 
the reference distribution's quadrupole-monopole tensor equal to the iden- 
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tity matrix q(y) = I. For a uniform density ellipsoid, the combined rotation 
R(t) and stretch A(t) transform it into the unit sphere, in units of i?o- 

The velocity vector for each particle in the inertial frame is U Q = (R t A+ 
R t A)y a . The kinetic energy of the rotating and vibrating system is 

T = E^ U «- U « 

a 

= (MRl/10) (-tv(A 2 n 2 ) + tr(i 2 )) . (46) 



The velocity V(t) is the tangent to the curve t i — > (R(t),A(t)) in the 
direct product space £0(3) x A 

V(t) = -3ln(t) + ]T Ofc(*)^-, (47) 

where A = diag(ai, a2, 03) and Q — RR 1 . The metric on the rotation group 
is extended to the direct product manifold so that d/dat is an orthonormal 
basis for the tangent spaces of A: 

g (R . A) (% h ^n 2 ) - -tr(ftiA 2 fi 2 ) (48) 
g(^)(^^n) = 0. 

Hence the kinetic energy Eq. (46) is proportional to the squared length of 
the velocity 

T = (MRl/10) S m) , A{t)) (V(t), V(t)). (49) 

The Lagrangian is the difference between the kinetic and potential en- 
ergies. The potential energy V(A) is a pure function of the axis lengths. 
The Euler-Lagrange equations may be obtained using a similar line of rea- 
soning as in the rigid body problem by enlarging the configuration space 
from 12 to 15 dimensions to accommodate the axis lengths. The Lagrange 
equations corresponding to the axis length vibrations are 

MRl „ _dJ^ _dV_ 
5 1 da* da>; ' 



The Euler equation Eq.(42) for rigid body motion is unaltered. 
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This dynamical system may be transferred from the direct product space 
SO (3) x A onto the ellipsoidal space Q. Consider the mapping 

(j) : SO{3) x A — > Q 

(R, A) i — ► q = <f>(R,A) = R l A 2 R. (51) 

Since any positive-definite symmetric matrix may be diagonalized by an 
orthogonal matrix and its eigenvalues are real positive numbers, the map- 
ping 4> is onto Q. But <f> is many-to-one for two reasons. First, although 
the positive eigenvalues of q G Q are unique, their ordering is not. For 
a fixed A g A, let K denote the group of all k g 50(3) such that 
kAk 1 g A. Then the transformation A i — ► kAk 1 permutes the entries 
of A = diag(ai, a 2 , a 3 ). If the entries of A are distinct, then K is the 
discrete subgroup of SO (3) with 24 elements consisting of the unit deter- 
minant matrices with exactly one nonzero entry, equal to ±1, in each row 
and in each column. Second, although the cigcnspaces of q are unique, the 
eigenvectors (the rows of R) are not. In the generic case of distinct eigen- 
values the rows are determined up to a factor of ±1 and permutations; left 
multiplication of R by k £ K corresponds to this indeterminacy. Define an 
action of the group K on the product space by (R, A) i — ► (kR, kAk*) for 
k g K. It is evident that q = <j>{R, A) = <f>(kR, kAk r ) for all k g K. 

The metric at (R, A) is invariant under the transformation group K. 
With respect to the if -action, the vector fields d/dai are permuted among 
themselves, and the right invariant vector field is transformed to 3i A d fc fi- 
Hence the metric of Eq.(48) is K- invariant. The AT-invariance of the metric 
implies that the Riemannian structure is well-defined on the manifold Q. 
The Lagrangian is also well-defined on the ellipsoidal space because the 
potential energy function is assumed invariant under permutations of its 
arguments, V(A) = V(kAk t ) for k E K. 

The model degrees of freedom corresponding to rotations and vibrations 
are coupled to each other in a simple way. As the body axis lengths vary, 
the moments of inertia change accordingly and, thereby, the rotation of the 
body is modified. But, at each instant in time, the moment of inertia is 
the rigid body value. For fluids, say a water droplet, the description must 
be enhanced so that the moment of inertia can assume the irrotational 
flow value. Physically this requires the inclusion of additional degrees of 
freedom corresponding to the internal vortex motion of the fluid. The 
Riemann ellipsoid model achieves that goal. 

3. RIEMANN ELLIPSOIDS 

The nine-dimensional configuration space of the Riemann ellipsoid model 
is the connected component of the general linear group, V = GL + (3, R). 
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The space V will be shown to be a principal bundle over the base manifold 
Q with structure group G. Locally the bundle V is the Cartesian product 
of Q and G, and points £ in V are defined by unique pairs £ = (q; S) for 
q € Q, S 6 G. As with the rigid rotor problem, the use of moving frames 
simplifies the geometrical analysis by allowing for tangent space bases that 
include the angular momentum and Kelvin circulation vectors. V is also 
a Ricmannian manifold whose metric is inherited from the kinetic energy 
on Euclidean space. The conservation laws for the angular momentum and 
Kelvin circulation are consequences of the metric's invariance with respect 
to left multiplication by the rotation group SO (3) and right multiplication 
by the structure group G, respectively. The dynamical equations for Rie- 
mann ellipsoids are derived from a Lagrangian that respects the metric's 
invariance properties. 

3.1. Bundle Structure 

The mapping n : V — ► Q , 7r(£) = is onto the space of all ellipsoidal 
shapes because cj> from SO (3) x A is onto Q and q — <j){R, A) — R t A 2 R = 

t(0 for Z = &A€V. 

7r is not one-to-one since V is nine dimensional and the ellipsoidal space 
Q is six dimensional. If both £i and £2 define the same ellipsoid, then 
7r(^i ) = 7r(£ 2 ) and (^ 1 <fi)(^ 1 ^i) t = /, or £1 i s an orthogonal matrix. 
Thus £1 and £2 define the same ellipsoid if and only if £1 and £2 are in the 
same left coset of G = SO(3), i.e., £1 e &G. There is a one-to-one and 
onto identification of the space of ellipsoidal surfaces with the left coset 
space 

Q = V/G = GL+ (3, R)/SO(3) . (52) 

By choosing a smooth set of left coset representatives, the space V may 
be identified locally with the Cartesian product of Q and G. Such a set is 
provided by a smooth mapping r from an open neighborhood U T C Q into 
V such that r(qi)G = r(q2)G if and only if q\ = (72, for qi, <?2 £W T . Given 
a smooth coset representative mapping, a diffeomorphism between U T x G 
and an open neighborhood V T of V is defined by 

U r xGcQxG — > V T CV 

(q;S) — > C = r(q)S (53) 

for q £ 1A T and S G G. The identification is not canonical since it de- 
pends on the arbitrary choice of the coset representative map r. The 
smooth choice of coset representatives r(q) = R l A for q = R t A 2 R is a 
well-defined one-to-one mapping only locally because the diagonalization 
of q is not unique. There is no smooth mapping r defined globally on Q. 
The space V is only locally diffeomorphic to the Cartesian product of Q and 
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G. Throughout this paper, whenever a chart for V is needed, a smooth 
choice of coset representatives of the form r(q) = R l A is taken in some 
open neighborhood IA T of Q, and the points of the bundle are given in an 
open neighborhood V T of V by £ = R t AS. This choice has the advantage 
of distinguishing among rotational, vibrational, and vortex motions. 

Consider the right multiplication of G on V: R g {£) = £<7 -1 for £ £ V 
and g £ G. G must be distinguished from the rotation group that acts 
by left multiplication on V . The point q £ Q corresponding to both £ 
and R g (C) is the same. Thus the transformation R g does not change the 
ellipsoidal boundary of the fluid, but it does change its internal motion. 
In terms of the local chart of Eq.(53) for V, the right multiplication is 
R g {q;S) = faSg' 1 ). 

The space V is a principal fiber bundle over the base manifold Q with 
structure group G, since the manifold V = V(Q, G, it) satisfies the following 
defining properties of a principal G-bundle [7] : 

• The space V is locally diffeomorphic to the Cartesian product of the 
base manifold and the structure group. Given a coset representative map, 
a bundle point in the corresponding local chart is £ = (q; S), where q £ Q 
and S £ G. Such charts are called local trivializations. 

• The projection n : V — ► Q is defined from the bundle onto the base 
manifold: 7r(£) = = q £ Q for £ = (g; S) £ V. The inverse image of 
a fixed point q £ Q in the base manifold is diffeomorphic to the structure 
group, TT^ 1 (q) = G, which is called the fiber over q. 

• A right multiplication R g : V — ► V is given by R g (£) — £g _1 and it 
satisfies n o R g = n. In terms of the local chart, R g (q; S) = (q; Sg~ x ) for 
t = (q;S)eP. 

3.2. Tangent Vector 

Consider a curve t i — > £(t) in the bundle V . Such a curve may be 
identified with the collective motion of a many-body system for which the 
trajectory of each particle is constrained by X Q (i) = £(i)y a , where y a 
is independent of time. The reference distribution y a is chosen so that 
its quadrupole-monopole tensor is the identity matrix, q(y) = I. With 
this choice the instantaneous dimcnsionlcss quadrupole-monopole tensor 
simplifies to q(X) = Hence the base manifold point 7r(£) is interpreted 
physically as the quadrupole-monopole tensor of the constrained many- 
body system. 

The velocity vector for each particle is U a = £y Q = ££ -1 X Q . Note that 
the velocity of each particle is a linear function of its position vector. The 
velocity vector may be expressed as the value of a right invariant vector 
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field on the group V at the point £, 

d 



V(t) = £4 



= -(DO C) for« = ^r 1 - (54) 

The vector field R u is right invariant with respect to the entire group V, 
in contrast to the invariant vector fields introduced in §2 for the rotation 
group. Although nominally derived for the rotation group, the results of 
§2 hold for any Lie group. This includes the concepts of left and right 
invariant vector fields and their transformation properties under left and 
right translation. The velocity u is an element of the Lie algebra Ms(R) of 
GL+(3,M), i.e., a real 3x3 matrix. 

A curve t i — ► (R(i) , A(t) , S (t)) in the direct product manifold 50(3) x 
AxG defines a curve t i — > f (i) = R{tf A(t) S \t) in V. The curve 1 1 — > S(t) 
in the structure group G defines the vorticity A = SS' 1 . A is an element 
of the Lie algebra g of the group G and must be an antisymmetric matrix. 
The curve 1 1 — > q(t) = 7r(£(i)) = R(tf A(t) 2 R(t) is in the base manifold Q 
and describes the rotation and vibration of an ellipsoid. 

However, any segment of a smooth curve t i — > ^(f) in the bundle V 
cannot be lifted uniquely to a smooth curve in the direct product manifold. 
Recall that the ordering of the elements of A is not unique. As discussed 
in §2, the ambiguity in the lifting is described by the subgroup K of the 
rotation group: for every k G K , £ = R t AS is invariant with respect to the 
transformation (R, A, S) i — ► (kR, kAk*, kS) in the direct product space. 

The velocity V(t) may be simplified to the sum of terms corresponding 
to rotational, vibrational, and vortex motions: 

u = iC 1 = R* (-V + A^A + AAA- 1 ^ R. (55) 

Each of these three vector fields on the bundle V may be expressed in the 
S0{3) xAxG chart (£ = R*AS) as 

(^RHm)i = -(%-ih = (VRh (56a) 



dRij / R 



(56b) 



A 



(%RtAAA-iRh = (^A)S - (AS)y ( — ) , (56c) 



-'y / s 
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and the velocity simplifies to 



V(t) = -Ql a ) R + -(Ba)s. (57) 



den 



A 



Note that (3?n)_R is a right invariant vector field on SO(2>) and (3?a)s is 
a right invariant vector field on G. This common notation for all right 
invariant vector fields should not be confusing since f2 is reserved for the 
Lie algebra S0(3) of SO(3) while A is reserved for the Lie algebra g of 
the vortex group G. A subscript £ indicates an invariant vector field on 
V evaluated at £. The subscript R (respectively, S) indicates an invariant 
vector field on SO (3) evaluated at R (respectively, on G evaluated at S). 

To prove Eq.(57), consider each of the three terms of u in Eq.(55) sepa- 
rately. For the rotational contribution to V(t), Eq.(56a), 

(%RHmh = ((R^RKhJ^ 
= (R'flAS)^^- 



(R*Q.) ik (AS)kj TTT— 
Ot;ij 

d 



dRki 

(ttfi)fl, (58) 



because, according to the chain rule, 



dRij dRij d£ki 

d 

= S m i5kj(AS) m i- 



d£ki 

Similar arguments apply to the vibrational and vortex terms, Eqs.(56b,56c) 



(AS)u^-. (59) 



3.3. Riemannian Structure 

The kinetic energy of the linear velocity field is 



a 

= (Mi2g/10)tr(&*). (60) 
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For X, Y in M 3 (R), define the metric at the point £ G V by 



g c ((^)i,(^) € )=tr(X(^*)y*). 



(61) 



This is a positive-definite bilinear form defined on each tangent space of V; 
hence, V is a Riemannian manifold. The kinetic energy is given in terms 
of the metric by 



Expanding the velocity into the three types of motion, Eq.(57), simplifies 
the kinetic energy to 



Note that the last term is a Coriolis coupling between the rotational and 
vortex degrees of freedom. 



Ricmann ellipsoid dynamics may be derived using the Lagrange formal- 
ism. The potential energy V is assumed to be a rotational scalar function 
on the base manifold, i.e., a smooth function V(A) of the axis lengths that 
is invariant under K, V(kAk t ) = V(A) for all k G K. The rotation, vibra- 
tion, and vortex motions of a Riemann ellipsoid correspond to the path in 
V that minimizes the action with respect to all smooth curves connecting 
two fixed points. 

The Lagrangian is invariant with respect to left multiplication by the 
rotation group 5*0(3) and with respect to right multiplication by the struc- 
ture group G. The potential energy is obviously invariant with respect to 
both of these multiplications. The kinetic energy is invariant because the 
metric is invariant: if r G 5*0(3), then left multiplication of the metric 
yields 



T=(MR 2 /10)g^(V(t),V(t)). 



(62) 




3.4. Lagrangian 



= tr ((rXr-^r^YirYr- 1 ) 1 ) 



(64) 



If g G G, then right multiplication of the metric gives 



= tr (Xttg-^g-yYt) 



(65) 
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By Noether's theorem, the invariances imply conservation laws. Left in- 
variance implies conservation of the angular momentum vector L, while 
right invariance requires conservation of the Kelvin circulation vector C. 
The Lagrangian is also invariant under time translation which implies con- 
servation of energy E = T + V . 

As with the rigid body problem, it is preferable for computational pur- 
poses to enlarge the configuration space V to the 27-dimensional Euclidean 
space whose coordinates are the 18 matrix entries of R and S, the six an- 
gular and vortex velocity components of Q and A, and the three entries 
of the diagonal matrix A [18]. This extension avoids the introduction of 
Eulcr angles for SO(3) and G, Moreover, two conservation laws follow be- 
cause the Lagrangian does not depend explicitly on the angular and vortex 
velocities. Since the allowed paths in the enlarged space must correspond 
to paths in the original configuration space V, they must satisfy eighteen 
constraint equations. In addition to the 5*0(3) constraints of Eq.(33), there 
are also constraints for the structure group: 

= S'S-I 

= SS* - A. (66) 

Eighteen Lagrange multipliers must be introduced in the form of real 
3x3 symmetric matrices a and r, and antisymmetric matrices \i and 
v. The Lagrangian L = L[R, R, Q, fl, A, A, S, S, A, A] on the twenty-seven 
dimensional configuration space is given by 

L = T(fl, A, A, A) — V(A) 

+\ tr[<T • (R l R - I)} + ±tr[r • (5*S* - /)] 

+§tr[/i • [RR 1 - 0)] + \tv[u ■ [SS 1 - A)]. (67) 

Except for the potential energy, the Lagrangian is just a quadratic function 
of the 27 coordinates and their derivatives. 

The conservation laws are derived by noting that L is independent of tl 
and A, and, therefore, the Euler-Lagrange equations for the angular and 
vortex velocities simplify to 

dT n dT 
fHj = -£ij/to — and v %j = - e ijk-^r- (68) 

The derivatives of the kinetic energy with respect to the angular velocity 
and vortex velocity are the vectors of angular momentum and circulation, 
respectively: 



L k = — = (MRl/h) [(of + a))u k - 2a iaj \k] (69a) 
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dT 

C k = -g^- = (MRl/h) [2a t a,Lu k - {of + a))\ k ], (69b) 



where i,j, k are cyclic. Note that the angular momentum depends on both 
the angular and vortex velocities. Because the two vectors are represented 
naturally by antisymmetric matrices, L — J^^k^k and C — J2Ck*k, La- 
grange's equations reduce to 

H = —L and v = C. (70) 

Lagrange's equations for the Rij and Sij are derived using arguments sim- 
ilar to Eqs.(38 - 40) and give the matrix equations 

fi=[Sl,fi] and v = [A, v] . (71) 

Hence, the dynamical equations for the rotational and vortex motion of a 
Riemann ellipsoid body are a Lax system 

^L=[Q,L] and jC=[K,C]- (72) 

Expressed in vector form, the precessions of the angular momentum vector 
and the Kelvin circulation vector are 

^ = —u) x L and ^ = -AxC. (73) 
dt dt y ' 

The two vector conservation laws arc 

£ (i2*L) =0 and ^ (S*C) = 0. (74) 

(Jib iXi-i 

For the axis lengths vibrations, Lagrange's equations in the variables a k 
yield the same formal equation (50), but the Riemann kinetic energy is 
substituted for the rigid rotor expression. The entire dynamical system for 
Riemann ellipsoids may be expressed as a finite-dimensional Lax system 
[19]. 



4. CONNECTIONS ON THE RIEMANN ELLIPSOID 

BUNDLE 

For many interacting systems, there are constraint forces in addition to 
those described by the potential energy V(A). The simplest case is the 
rigid body for which the vortex velocity vanishes A — 0. This holonomic 
constraint simplifies the dynamical system to the rigid body theory of §2. 
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But constraints are not typically holonomic. For example, an irrotational 
fluid (like a water droplet) has zero circulation, C = 0. Another example 
is the falling cat [20], for which the angular momentum vanishes, L = 0. 
In these cases the vortex velocity is proportional to the angular velocity [3] 

irrotational flow: Afe = ^ fl » a j (75a) 

aj + aj 

a i + a i 

falling cat: Afe = — -u>k, (75b) 

where i, j, k are cyclic. For a holonomic constraint the structure group dy- 
namical equation, 5 = AS, may be solved explicitly for S(t). For example, 
if A = J2^k?k is constant, then S(t) = exp(iA)5(0), where 5(0) denotes 
the value of the structure group element at t — 0. Unless the axis lengths 
and angular velocity vector are constant in time, the irrotational flow and 
falling cat problems arc not generally intcgrable systems. The typical non- 
holonomic constraint for a Ricmann ellipsoid is a proportionality between 
the vortex and angular velocity components, Afe = f^{A)u>^ with a factor 
fk(A) that depends on the axis lengths. 

Suppose a nonholonomic constraint determines the vortex velocity in 
terms of the axis lengths and the angular velocity vector. A curve 7 : 1 1 — ► 
q(t) in the base manifold Q defines curves 7 : t 1 — ► £(t) = (q(t);S(t)) 
in the bundle V known as horizontal lifts. To obtain the horizontal lift, 
take the unique tangent vector field, Eq.(47), of the curve 7 in the base 
manifold, and express the vortex velocity K(t) in terms of the angular 
velocity Q(t) and vibrational velocity A(t) by applying the nonholonomic 
constraint. Eq.(57) defines a vector field on the bundle that integrates 
to the horizontal lift 7 through the initial point 5(0) 6 G. The bundle 
vector field is only defined on the subbundlc n^ 1 (q(t)) lying above the base 
manifold curve. Note that the horizontal lift depends on the association of 
a unique vector field on the bundle to a vector field on the base manifold. 
As it will be seen later, this association is the differential structure, or 
connection, on a bundle. 

The concept of a horizontal lift is physically natural. It says that a many- 
body system responds to rotations and vibrations (described by a curve 7 
in the base manifold) by internal vortex motions (described by a curve 7 
in the bundle). This response is determined typically by a nonholonomic 
constraint that depends ultimately on the nature of the forces between the 
particles. The mathematical description of the response is a horizontal lift 
7 in the bundle, see Figure 3. 

The horizontal lift is determined by a connection in the principal bundle. 
Since the theory of connections is the foundation of differential geometry, a 
brief review will be given here. The theory of Riemann ellipsoids is regarded 
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as a paradigm for the general gauge theory of two interacting sets of physical 
degrees of freedom, e.g., quadrupole and pairing [21], symplectic[22] and 
interacting boson[23], etc. Although the theory of bundle connections is 
well known in particle physics, it is less familiar in many-body science. 

Suppose V = V(Q, G, n) is a principal bundle over the base manifold Q 
with structure group G and right multiplication R g (£,) = £ • g for all £ G V 
and g G G. R g is assumed to be a group homomorphism, R gi R g2 — R gig2 , 
or (£ ' gi) ' .92 = £ • (.91.92)- The Riemann ellipsoid bundle V = GL + (3,M) 
is an example, but the theory of connections will be cast in its general 
setting. In a local trivialization, a bundle point is £ = (q; S), where q G Q 
and S € G, the projection map is ir(q; S) = q, and the right multiplication 
is R g (q;S) = (q;Sg- 1 ). 

For each point £ in the bundle, denote the tangent space to V by T(P . 
The vertical space is a subspace of T{P consisting of the tangents to 
curves t 1 — ► (q;S(t)) in the fiber (q G Q is fixed). It may be defined 
without reference to a local trivialization by 

V 6 = {X G T(P I ir*X = 0} . (76) 

Among the vertical vectors, the left invariant vector fields on the struc- 
ture group are singled out for special consideration. Suppose that the 
structure group is a Lie subgroup of a real matrix group GL(d, M), and 
the Lie algebra g of the group G may be identified with an algebra of 
matrices. If A S g is a Lie algebra element, then the fundamental vec- 
tor field, denoted by A*, is the left invariant vector field La on the fiber 
G. If {f , a = 1, . . . ,p} is a basis for the p-dimensional Lie algebra, then 
{(f„)s, a = 1, . . . ,p} is a basis for the tangent space to the structure group 
at the point S G G. The fundamental vector fields define a Lie algebra ho- 
momorphism, [Ai, A 2 ]* = [A*, Aj], for Ai, A 2 G g. The right multiplication 
in the bundle satisfies tt o R g = tt and, hence, 7r* o (R g )* = 7r*. Thus (R g )* 
is a linear transformation from the vertical space onto the vertical space 
V^. g , and, according to Eq.(26a), it is given by (R g )*(A*)s = (Ad g A)g g _ 1 . 

A connection is a smooth assignment F of a subspace of the tangent 
space T(V to each point £ G V. is called the horizontal subspace. A 
vector X G T(P is called vertical, respectively horizontal, if it lies in Vjt, 
respectively in H^. To be a connection, the assignment must satisfy two 
properties: 



(1) Tf.T = H t V ( (direct sum) 

( 2 ) H H-g = i R g)* H £ ■ 



(77a) 
(77b) 
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The first condition means that every tangent vector X E T(V can be 
uniquely written as 

X = Y + Z, where Y E and Z E V ( , (78) 

and Y, respectively Z, are called the vertical, respectively horizontal, com- 
ponents of X. The horizontal projection Y E of X E T(V is denoted by 
Y = hX. 

Since the kernel of 7r» at £ E V is the vertical subspace V^, the image 
of 7r» is T q Q, and the tangent space is a direct sum of the vertical and 
horizontal subspaces, the linear transformation 7r» is an isomorphism from 
the horizontal subspace onto the tangent space of the base manifold 7r» : 
— > T q Q , where <z = tt(£). If T E T q Q is a tangent vector to the base 
manifold, then its horizontal lift is the unique horizontal vector T E such 
that 7r*T = T. Given a basis of smooth vector fields to the n-dimensional 
base manifold at q E Q, {(ei) q ,i — 1, . . . ,n}, their unique horizontal lifts 
are denoted by {(E,)|, i = 1, . . . , n}, where 7r(£) = q. The set {(Ej)j} is a 
basis for the horizontal subspace and 

(E i )| = (e i ) g -^r«(0(O S , (79) 

a 

where, in a local trivialization, £ = (q;S). The coefficients Tf must be 
smooth real- valued functions on the bundle V '. 

The second defining property of a connection, Eq.(77b), asserts that 
(Ei)jr.g = (i?g)„(Ei)^. In particular, when £ = (q; e), where e is the group 
identity and g = S^ 1 E G, the right translation of a horizontal basis vector 
at the structure group identity, (Ej)( 9;S ) = (i? s -i)*(Ej)( ?;e ), is 

(Ei)(, iS ) - (e^-^rjfoeXAds-ifJS 

a 

= (e^+^Tr^X^Js, (80) 

a 

where Eqs.(19, 26a) are used. Hence, the right invariance of the horizontal 
spaces implies that the functions T arc determined by their values at £ = 
(q; e). The functions T?(q) = Tf(q; e) are called Christoffel symbols. 

Suppose q(t) is a smooth curve in the base manifold Q and its tangent 
is given by T(t) — ^(*)( e »)g(t)- The tangent to the horizontally lifted 
curve 1 1 — > (q(t);S(t)) is f — ^2T % (t)CEi)^ t y iS (t)), and it may be decom- 
posed into horizontal and vertical components 



(81) 
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where Sjk are the matrix entries in GL(d,M). The components Sjk are 
given by 



S jk = T(Sjk) 



i a 



-E r wE r "(^))(f a ^, 



(82) 



a 



where Eq.(16b) is applied. In matrix form the differential equation for the 
horizontally lifted curve is 



If the term in parentheses is independent of time, then the differential 
equation integrates to an exponential. 

This finishes the required mathematical preliminaries. The theory will be 
applied now to the Riemann ellipsoid model. Let ti denote, as before, the 
natural basis for the algebra of 3 x 3 antisymmetric matrices, {ti)jk = £ijk- 
A basis {(ei) q ,i = 1,...,6} for the tangent space at q e Q consists of 
the three right invariant vector fields (ej) g = (ftejfl on the rotation group 
S*0(3) and the three vibrational vector fields (ei + 3) 9 = (d/dai)A on A for 
i = 1,2,3. A tangent vector to a curve in the base manifold is given by 
Eq.(47) or 



A basis for the Lie algebra Q of the structure group is the set {f^ = C,, i = 
1, 2, 3}. The horizontal lift of the curve is given by Eq.(57) or 

3 




(83) 



3 




(84) 




8 = 1 




(85) 



Therefore the tangent vectors 




(86a) 



d_ 



(86b) 
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for i = 1,2,3, are a basis for the horizontal subspace. They show that 
the Christoffel symbols for Ricmann ellipsoids vanish for the vibrational 
vectors and simplify to a diagonal form for the rotational vectors 

r?( ff ) = S t g) . (87) 

In particular, the special rotational modes correspond to the following 
Christoffel symbols: 

{0 rigid 
2ai<ij / (a? + a?) irrotational (88) 
(aj + a'j)/(2a i aj) falling cat, 

where k are cyclic. Note that the irrotational Christoffel symbols are 
at most one and the falling cat's are at least one. 

4.1. Riemannian Connection 

The irrotational flow and falling cat connections arise naturally from the 
Riemannian structure and the group structure, respectively. The horizontal 
subspace H$ for irrotational flow is defined as the orthogonal complement 
to the vertical subspace . Denote the vector space of all 3 x 3 symmetric 
matrices by m. The orthogonal complement is given explicitly by 

H 6 = {(£y)e £ T(P rem}. (89) 

To prove this, suppose (&a)s, A g g, is a vertical vector and (31y)z, Y S m, 
is a horizontal vector. These two vectors are orthogonal since, by Eq.(56c), 

g c ((3U) S ,(3V)«) = gt((%IPAAA-iRh,(%Yh) 

= ti(& AAA' 1 RitftY*) 
= tx(R t AAARY) 
= -tr (Y t R t AAAR) 

= 0. (90) 

In addition, observe that the sums of the dimensions of the vertical space 
and the horizontal space add to the dimension of the tangent space T(P . 
Therefore the tangent space is a direct sum of the horizontal and vertical 
subspaces. Moreover, if (9fy){ is a horizontal vector and g £ G, then right 
invariance implies 



(i? g )*(fty)£ = (3ly)| g -i, 



(91) 
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or (Rg)^H^ = H^ g -i. Since the assignment of the horizontal subspace 
is also smooth, it defines a connection L on V . The Christoffel symbols are 
found by determining explicitly the horizontal lifts of the vibrational and 
rotational vectors. 

The vibrational vectors are horizontal since Y = R t A~ 1 AR is a sym- 
metric matrix, sec Eq.(56b). Therefore the Christoffel symbols for the 
vibrational vectors are zero. But the rotational vectors arc not horizontal 
because 

8e((ttc,)R, (ttejs) = tr (At t Ae b ) = -25 ib aja k ^ (92) 

for i, j, k cyclic. Note that the inner product of two vertical vectors is also 
nonzero, 

g^ejs, (B c Js) = ~ tr ( A ^at b ) = + 4) (93) 

for a,j, k cyclic. In order for (Ej)^ to be the horizontal lift of (R^Jr, it is 
necessary and sufficient that, for b = 1, 2, 3 , 

= ^((EOf.CRcJs) 

= g€((K ei )B + r?(g)(3le a )s,(KeJs) 

= -2^a, afe +r 4 b (g)(a2+ a 2). (94) 

The off-diagonal Christoffels for the rotational vectors vanish, and the 
diagonal values are 

rj(g) = JT*^ cyclic). (95) 

Thus, the Riemannian connection for which the horizontal space is perpen- 
dicular to the vertical space corresponds to irrotational flow. 

4.2. Invariant Connection 

The falling cat connection is the invariant connection on the Lie group 
V. Since Q is the algebra of antisymmetric matrices and tn is the vector 
space of symmetric matrices, the Lie algebra of the group V is a direct 
sum of vector spaces, M 3 (R) = Q © m. Moreover the vector space m is 
invariant with respect to the adjoint group transformation, Ad g (m) C tn 
for all jgG. These two properties of tn are necessary and sufficient for 

Hi = {(£y)c = -(%AdtY)t G T ( V I Y g m} (96) 

to be a horizontal subspace [7]. In order to see that, note that the vertical 
vectors can be expressed in left invariant form, 



V 6 = {(X A )s = 



-(£s«as) € e T 6 V I Aeg}. 



(97) 
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The tangent space to the bundle at £ is a direct sum of the horizontal 
and vertical subspaces, because every matrix is a linear combination of 
a symmetric matrix Y and an antisymmetric matrix S AS. The right 
invariance of the horizontal subspaces is a consequence of 

{R g )*{$.Ad ( Y)t = = (^Ad £9 _ 1 .4d 9 y)c 9 -i e H ig -i, (98) 

since Ad g Y € m for all g € G and Y <E m. The assignment of the subspaces 
is smooth, so is indeed a horizontal subspace. 
The relation 



(99) 



shows that the vibrational vectors are horizontal (S t A 1 AS is symmetric), 
but the rotational vectors are not horizontal since 

(#e,k = -(^Ji*Ciij)« = {&stA-HiAs)i (100) 

and S t A~ 1 z i AS is not symmetric. If the matrix S t A^ 1 t i AS is expressed as 
a sum of symmetric X s and antisymmetric X a parts, i.e., X s — [S t A~ 1 tiAS— 
S t At l A~ 1 S)/2, X a = {S'A^t.AS + S t At l A- 1 S)/2 1 the angular momen- 
tum may be written as a sum of horizontal and vertical vectors 

(ttejii = {&x.)i + {^x a )i eH^O V^. (101) 
The horizontal lifts of the angular momentum vectors are 

(E^ = (X tt ) R + Tt(q)(3lt a )s 

= (Hx.h + ((fix a )t-mq)(/is*t a s)t), (102) 

where (£>x s )$ is the horizontal lift and the two vertical vectors in the 
parentheses must cancel. Therefore, the invariant connection is given by 

TUq)*a = (A-hiA + At l A- 1 )/2, (103) 

or the Christoffel symbols are diagonal and 

TM = <|±^, (104) 

where i,j, k are cyclic. 

4.3. Curvature 

Given a bundle connection T, define a 1-form u> on V with values in the 
Lie algebra Q of the structure group G as follows: if X G T^P is a tangent 
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vector in the bundle and Z is its vertical component, then there is a unique 
element A e such that the vertical component is the fundamental vector 
field A* = Z. At i e V, define u{X) = A. The form u is called the 
connection form of the connection T. The connection form satisfies three 
general properties: 

(1) u>(X) = 0, if and only if X is horizontal; 

(2) w(A*) = A for every Aeg; 

(3) u)((R g )*X) — giv(X)g~ 1 for every g £ G and every vector field X. 

The first two conditions arc an immediate consequence of the definition 
of the connection 1-form. The third property is proven by considering 
separately horizontal and vertical vectors. If X is horizontal, then (R g )*X 
is also horizontal and both sides in property three are zero. If X is a vertical 
vector at £ G V, then it equals some fundamental vector field, X^ = (A*)^ 
and the third property is a consequence of Eq.(26a). Conversely, if the 
second and third conditions are satisfied by some 1-form u), then it defines 
a connection for which the horizontal subspaces are the vectors X satisfying 



If u> is the connection 1-form, then its exterior derivative du) is a 2-form, 
i.e., a bilinear, antisymmetric form defined on pairs of tangent vectors X, Y 
to the bundle by 



where hX and hY are the horizontal components of the vectors X and Y, 
respectively. If either X or Y is vertical, then £l(X, Y) = 0. If both X and 
Y arc horizontal, then 2Q.{X,Y) = -u([X,Y]). 

The curvature form of the Riemann ellipsoid bundle is, for 1 < i, j < 3, 
as follows: 



u>(X) = 0. 



dw{X, Y) - \ {X(co(Y)) - Y(oj(X)) - u([X, Y})} . 



(105) 



The curvature form Q of uj is defined by Q — {du)h, 



n(X,Y)=duj(hX, HY), 



(106) 




dai ' ddj 







(107a) 




%j (Ads-itj) 
e tj k9\k {Ad s -itk) 



(107b) 



(107c) 



where the field tensors are given in terms of the connection by 



(108a) 



2 dai 
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^-rjrj) 



(i,j,k cyclic). 



(108b) 



To derive the curvature expression, use the commutation relations among 
the horizontal vectors: 



d d 
'aj ' daj 

d 

da~^ 







dT 
da 



(109a) 



(109b) 



[e,,e,] = - SiJk (E k + (r k k -r\r>)(Ad s -itkr). (i09c) 



In general, if is any g-valucd p-form, then its covariant derivative D9 
is a 0-valued (p + l)-form given by the composition of horizontal projection 
and exterior differentiation, D8 = (dff)h. The covariant derivative DO and 
the exterior derivative d6 are totally antisymmetric in their arguments that 
are p+ 1 tangent vectors to the bundle. The curvature form is the covariant 
derivative of the connection 1-form, il = Dm. The Bianchi identity states 
that the covariant derivative of the curvature form is zero, DQ = 0. For any 
three horizontal vectors X, Y, Z, the covariant derivative of the curvature 
form is given by 

3Dn(x,Y,z) = x(n(Y,z)) + Y(ri(z,x)) + z(n(x,Y)) 

-n([x, y],z)- n([Y, z],x)- Q([z, x],y). (no) 

If all three horizontal vectors are vibrational, then it is clear that the covari- 
ant derivative is zero. Although not self-evident, the covariant derivative 
is also trivially zero if the three horizontal vectors are rotational. If two 
arguments are vibrational and one rotational, then 

"^^"w'-Kl 1 -©^-" 1 <»•) 

for all 1 < i, j, k < 3. If one argument is vibrational and the other two 
rotational, then 

d ( dH \ 

= 3 £>fi(Ei, Ej, — ) = e ijm f rjT fcj + rjXw - % km + -g^j (Ads-itm). 

(112) 

The latter equation is proven using the relations 



{$. ti )s(Ad s -iej) = Ad s -i[ti,tj\. 



(113a) 
(113b) 
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Conversely, suppose the Bianchi identities, Eqs.(lll,112), are satisfied by 
the field tensors. Define three "displacement vectors" , = (%u, ^2i, ^3i), 
i.e., the i th column of the tensor is the vector D^. The first Bianchi identity, 
Eq.(lll), is equivalent to the vanishing of the curl of each displacement 
vector 

V„xD, = 0. (114) 

Since A is a simply connected region, each irrotational displacement vector 
is the gradient of some "electric potential" , 

Di = ^V o 0<. (115) 

By substituting the potentials fa for X, the second Bianchi identity, Eq.(112), 
simplifies to 

= e ijm { A (2 * m _ ^ + Mj) + {Tl _ ^ + {T i _ ^} 

x(Ads-itm)- (116) 

If the potentials arc set equal to the Christoffel symbols, fa =T\, then the 
above equation determines, up to additive constants, the 91 tensor. 

4.4. Projection to the Base Manifold 

The connection one-form ui and the curvature two-form Q on the bundle 
V can be projected down to the base manifold. The projection depends 
on the choice of a local trivialization r : Q — ► V . If 9 is any g-valued p- 
form on the bundle, then its pullback t*8 is a g- valued p-form on the base 
manifold. The pullback of the connection one-form A = t*oj is a one-form 
on Q given by 

A q (X) = u(t*X) for any tangent vector X eT q (Q). (117) 

For our usual choices of a local trivialization, r(q) — and an adapted 
basis of tangent vectors, {(ej) g , 1 < i < 6}, the pullback is just 

A(-) - (118a) 

A(R tk ) - T k k {q)t k . (118b) 

If the pullback is applied to the tangent vector of a curve in the base 
manifold, then it produces the vortex velocity 

A(V(t)) - -uj k T k k t k = -A(i). (119) 
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The pullback of the curvature is defined similarly, F = r*f2, and it is a 
two-form on the base manifold: 

F q (X,Y) = Q(t*X,t*F) for all tangent vectors X, Y e T g (Q). (120) 

The tensor F is evaluated to be 

oai oaj 

n 'L i ^ ) = %jtj (i2ib) 

F(R ti ,R t .) = EijkMktk. (121c) 



The Bianchi identity on the bundle takes a slightly more complicated 
expression on the base manifold 

dF + [A, F] = 0, (122) 

where dF denotes the exterior derivative of F, see Eq.(llO), and the com- 
mutator is the three-form defined by 

3[A,F}(X,Y,Z) = [A(X),F(Y,Z)} + [A(Y),F(Z,X)] + [A(Z),F(X,Y)} 

(123) 

for all tangent vectors X,Y, Z e T ? (Q).The identity Eq.(122) is verified 
by letting X, Y, and Z run over the basis of the tangent space. 



5. CONCLUSIONS 

The connections corresponding to rigid rotation, irrotational flow, and 
the falling cat were shown to be natural geometrical or group-theoretical 
concepts. Although not mathematically natural, other choices of Christoffel 
symbols define nonholonomic constraint forces that are not excluded by 
physical law. For example, the S'-type Riemann ellipsoids are a sequence of 
special case solutions for which the angular momentum, Kelvin circulation, 
and the angular and vortex velocity vectors are aligned with a principal 
axis, say the 1-axis [2, 3]. This sequence is indexed by a continuous real 
parameter / restricted to the interval [—2, 0]. There is only one horizontal 
lift to consider and the Christoffel symbol is given by 



(124) 



34 



G. ROSENSTEEL AND J. TROUPE 



At / = 0, the connection yields rigid rotation, and, at / = —2, it is 
irrotational flow. The S-type ellipsoids are the simplest models that allow 
for a continuous interpolation between rigid rotation and irrotational flow. 
This connection has no natural geometrical or group-theoretic significance 
- but it does model a variety of rotating physical systems. 

The unity found in rotating systems is remarkable [24, 25]. Besides 
astrophysical systems [2, 3] , the Riemann ellipsoid model applies to classical 
liquid droplets [26] and quantum atomic nuclei [5, 27, 30, 31]. For example, 
the kinetic energy of a Riemann ellipsoid, Eq.(63), can be derived from self- 
consistency of a mean field approximation to the quantum fcrmion system 
of rotating nucleons [28] . 

A basic science problem is to determine the connection from the interac- 
tions among the particles that form a rotating system. The horizontal lift of 
a curve describing the rotation of a body depends on the coupling between 
the collective degrees of freedom of the Riemann ellipsoidal model and 
those of the individual particles. Equations that incorporate these interac- 
tions into the gauge theory and whose unique solution are the Christoffcl 
symbols are required for a complete theory of collective rotation. They 
must involve a coordinate independent object and the curvature form is 
the obvious candidate. The Bianchi identity, DQ = 0, partially determines 
the Christoffel symbols, but it is not sufficient. There must be another 
equation that relates the bundle curvature to the microscopic physics. 

The analogy with classical theory of electricity and magnetism is helpful. 
This is an Abclian gauge theory with space-time as the base manifold and 
U(l) as the structure group. The curvature form is the Faraday tensor and 
the Bianchi identity is equivalent to two of Maxwell's equations, V • B = 
and V x E + dH/dt = 0. The other two of Maxwell's equations involve 
the current and are equivalent to a second independent equation for the 
Faraday tensor. It is the analogous second equation for the gauge theory 
of rotating Riemann ellipsoids that is unknown. 

The quantization of a gauge theory is achieved by constructing associated 
bundles. For each irreducible representation of the structure group G ~ 
SO(3), an inequivalent associated bundle and a different quantization can 
be defined. In a subsequent article, the associated bundle theory will be 
presented. The quantized gauge theory applies directly to nuclear structure 
physics, e.g., the liquid drop model [29] is the associated bundle constructed 
from the one-dimensional representation of the structure group. Although 
the quantization of the Riemann model was considered previously [32, 33], 
the differential structure of the quantized theory has not been studied yet. 
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FIG. 2. Right translation (R g )* of a left invariant vector field £jq. 
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Principal G-Bundle 




FIG. 3. A curve 7 in the base manifold lifts to a curve 7 in the bundle. The curve 
in the bundle is called the horizontal lift. It depends uniquely on the bundle connection. 



